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Abstract 

We present work going on calculating the kaon Bk parameter in the 1/N C 
expansion. The goal of this work is to analyze analytically in the presence 
of chiral corrections this phenomenologically very important parameter. 
We present the method used and preliminary results for the chiral limit 
value for which we get B K = 0.29 ± 0.15. We also give some analytical 
indications of why the large N c prediction of Bk = 3/4 may have small 
1/N C corrections in the real case. 
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1 Introduction 



Indirect Kaon CP-violation in the Standard Model (SM) is proportional to the 
matrix element 



(K°\K°) = -iC AS= 2C(v)(K°\Jd 4 yQ A s =2 (y)\K c 



16 - 

-iC A s=2 y B K f K m 2 K 



with 



Qas= 2 (x) = AL^x)L^(x) ; 2L^x) = [s^(l - 75 )d] (x) . 



(1) 



(2) 



and C(v) is a Wilson coefficient which is known in perturbative QCD at next- 
to-leading (NLO) order in a = as /it in two schemes [1], namely, the 't Hooft- 
Veltman (HV) scheme (MS subtraction and non-anti-commuting 75 in D 7^ 4) 
and in the Naive Dimensional Regularization (NDR) scheme (MS subtraction 
and anti-commuting 75 in D 7^ 4). The coefficient C A s=2 collects well known 
functions of the integrated out heavy particle masses and Cabibbo-Kobayashi- 
Maskawa matrix elements. For comprehensive reviews where these factors and 
complete details can be found see [2]. The Wilson coefficient C(p) is 



C(u) = l + a(z/) 



72 
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(3) 



where 71 is the one-loop AS = 2 anomalous dimension 

7l = ^( 1 "^) ; 
72 is the two-loop AS = 2 anomalous dimension [1] 
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(5) 



and Pi and /3 2 are the first two coefficients of the QCD beta function. 

The so-called kaon parameter defined in (1) is an important input for 
the unitarity triangle analysis and its calculation has been addressed many times 
in the past. There have been four main techniques used to calculate the B K 
parameter: QCD-Hadronic Duality [3, 4], three-point function QCD Sum Rules 
[5], lattice QCD and the 1/N C (N c = number of colors) expansion. For a recent 
review on the unitarity triangle where the relevant references for the inputs can 
be found see [6]. For recent advances using lattice QCD see [7, 8]. 
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Here, we would like to present work going on determining the Bk parameter 
at NLO order in the 1/N C expansion. That the 1/N C expansion would be useful 
in this regard was first suggested by Bardeen, Buras and Gerard [9] and reviewed 
by Bardeen in [10, 11]. There one can find most of the references to previous 
work and applications of this non-perturbative technique. 

Recent related work to the one we will discuss in the following sections can 
be found in [12, 13] where a NLO in 1/N C calculation of Br within and outside 
the chiral limit is presented. There, the relevant spectral function is calculated 
using the ENJL model [14] at intermediate energies while at very low energies 
and at very large energies the chiral perturbation theory (CHPT) and operator 
product expansion OPE) results,respectively, are used. Another calculation of 
B K in the chiral limit at NLO in the 1/N C is in [22]. There, the relevant spectral 
function is saturated by the pion pole and the first rho meson resonance -minimal 
hadronic approximation (MHA). In [23], the same technique as in [22] was used 
but including also the effects of dimension eight operators in the OPE of the 
AS = 2 Green's function and adding the first scalar meson resonance to the 
relevant spectral function. 

2 Technique 

All the details on the A-boson method were given in [13]. In particular, in that 
reference it was shown how short-distance scale and scheme dependences can be 
taken into account analytically in the 1/N C expansion. Here, we just sketch the 
procedure introducing the notation. We want to calculate the two-point function 
[12, 13] 

n A 5= 2 (g 2 ) = i J d A e* x (0|T (P^(0)P K »( X ) e*™) |0> (6) 

in the presence of the long-distance AS* = 2 effective action of the Standard Model 
Ta5=2- After reducing the kaon two-quark densities the two-point function (6) 
provides the matrix element in (1). 

The effective action T A 5=2 reproduces the physics of the SM at low energies 
by the exchange of a colorless heavy AS = 2 A-boson. To obtain it [13], we make 
a short-distance matching analytically between 

Ta5=2 = -C AS=2 C(u) J d 4 yQ AS=2 (y)+h.c. (7) 

and -in our approach- 

Tld = 2 g A s=2^c, ■■■) j d A yX^y) L^y) + h.c. (8) 

Where, we have chosen to regulate the long-distance effective action in four di- 
mensions with a cut-off /ic- These choices are perfectly compatible with keeping 



2 



the short-distance scale and scheme independence analytically exact and we be- 
lieve that at low energies, where the relevant degrees of freedom are not quarks 
and gluons but hadronic degrees of freedom, are more natural. 

As said above, this matching takes into account exactly all the short-distance 
scale and scheme dependences as well as the choice of evanescent operators. We 
are left with the coupling of the X-boson long-distance effective action completely 
fixed in terms of the SM ones 



Ml 



1 + a (71 log {^^j + Ar 



(9) 



The one-loop finite term Ar is scheme dependent 



Ar 



NDR 



8 V NJ 8 V NJ 



(10) 



and makes the coupling |<7as=2| scheme independent to order a 2 . This coupling 
is as well analytically scale-z/ independent at the same order. Notice that there 
is no dependence on the cut-off scale jic -this feature is general of four-point 
functions which are product of conserved currents [13]. 

In the procedure described above, one also produces the standard leading and 
next-to-leading resummation to all orders of the large logs in [as log(Mw/v)] n 
and a s [as \og{M w /v)] n . And this has been done in the two-schemes described 
before; namely, NDR and HV. 

Once the long-distance effective action Tld is fully fixed, we are ready to 
calculate the relevant matrix element. 



(if°(g)|e irAS = 2 |K (g)> 



(ir( q )\^\K ( q )) 
d4 Px 9as=2 
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-iC AS = 2 Y BrQ f> 



(11) 



where q 2 is the external momentum carried by the kaons. The basic object is the 
four-point function 



n 



(p 2 x , q 2 )=t 2 4(K°( q )\ Jd 4 x j ^ye-^^T {L»{x)L»{y))\K\q)) . (12) 



At large N c , this four-point function factorizes into two two-point functions at 
all orders in quark masses and external momentum q 2 . The disconnected part of 
the four-point AS* = 2 function is 



g^n^ isconn (p^, q 2 ) 



(2n)^\p x )8q 2 l 



(13) 



which leads to the well-known large- N c prediction 



R Nc - - 
4 



(14) 
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At next-to-leading in the l/N c expansion, one has 

1 



B 



K 



poo 

/ dQ 2 F[Q' 
Jo 



with Q 2 the X-boson momentum in the Euclidean space and 



F[Q 2 ] = 



_ 9AS=2 

8tt 2 , 



lim 



Q 2 9llu U^ nn XQ 2 , q 2 ) 



Q 2 + M 2 X 



(15) 



(16) 



The next point is the calculation of (16). There are two energy regimes 
where we know how to calculate F{Q 2 ) within QCD. The first one, it is when 
Q 2 » 1 GeV 2 while q 2 is kept small since we will have to put it on the kaon 
mass-shell. In this regime, using the operator product expansion within QCD 
one gets 



^n^ n .(QV) = E 



C^,Q 2 )(jT(q)\Q^ +2 \K%q)) 



(0 



n=2 



Q 



2n 



(17) 



where Q,2n+2 are local AS = 2 operators of dimension In + 2. In particular, 



Qe 



4 J d 4 xL^{x)L fl {x) 



and 



with 



C 6 (v,Q 2 ) = -87r 2 7l a 



Fi 



1 + a 



(A - 71) log 



+ 0(a 2 ) 
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Ar 
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(18) 



(19) 



(20) 



where the a 2 term was not known before. The finite term Fi is order N c and 
therefore this a 2 term is of the same order as the leading term. In fact, at the 
same order in N c there is an infinite series in powers of a. 

For the list and a discussion of the dimension eight operators see [23, 25]. 
In [23] there is a calculation in the factorizable limit of the contribution of the 
dimension eight operators. Numerically, the finite term of order a 2 competes with 
that result when Q 2 is around (1 ~ 2) GeV 2 . 

The second energy regime where we can calculate -F[Q 2 ] model independently 
is for Q 2 — > where the effective quantum field theory of QCD is chiral pertur- 
bation theory. The result is known [13, 22] up to order p 4 



F*{Q 2 



4a (2£i + 5L 2 + L 3 + L 9 ) Q 2 + 



(21) 



with F the chiral limit of the pion decay constant f n = 92.4 MeV. 

We still need the intermediate energy region for which we use the large iV c 
hadronic model described in the next section. 
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2.1 Large N c Hadronic Model 

Up to this point, all the results are model independent. In particular, we have seen 
in the previous section that there are two energy regimes which can be calculated 
within QCD. In this section, we describe a large N c hadronic model which will 
provide the full YI^ nn (Q 2 ,q 2 ) which contains these two regimes analytically. 

The large N c hadronic model we use was introduced in [24]. The basic objects 
are vertex functions with one, two, . . . two-quark currents or density sources 
attached to them, referred to as one-point, two-point,. . . vertex functions. These 




Figure 1: One-point, two-points, three-points, • • • vertex functions. The crosses 
can be vector or axial-vector currents, scalar or pseudoscalar densities. 

vertex functions are glued into infinite geometrical series with couplings gy for 
vector or axial-vector sources and g$ for scalar or pseudo-scalar sources. In these 
way one can construct full n-point Green's functions in the presence of quark 
masses -see for instance, how to get full two-point functions in Figure 2. 




Figure 2: Infinite geometrical series which gives full two-point functions at large 
N c . The crosses that glue the vertex functions are either gy for vector or axial- 
vector sources or g s for scalar and pseudo-scalar sources. 

The basic vertex functions in Figure 1 have to be polynomials in momenta 
and quark masses to keep the large N c structure but just keeping the first octet of 
hadronic states per channel, i.e., the pseudo-scalar pseudo-Goldstone bosons, the 
first vectors, the first axial- vectors and the first scalars non-dynamically generated 
resonances. The coefficients of the vertex functions are free constants of order N c . 
Imposing chiral Ward identities on the full Green's functions one can fix most 
of these unknowns. Chiral perturbation theory at order p 4 and the operator 
product expansion in QCD help to fix the rest of the free coefficients of the 
vertex functions. 

The full two-point Green's functions obtained from the resummation in Fig- 
ure 2 agree with the ones of large N c when one limits the hadronic content to 
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be just one hadronic state per channel -our model does not produce any less or 
more constraints and all parameters can be fixed in terms of resonance masses 
[13]. Introducing two or more hadronic states per channel systematically could 
be done -it would just be much more cumbersome. To our knowledge, all the 
low-energy hadronic effective actions used for large N c phenomenology are in the 
approximation of keeping the first resonances below some hadronic scale and in 
some cases not in all channels. 

All the process of obtaining full Green's functions can be done in the presence 
of quark masses. In fact, in [13] two-point functions were calculated outside the 
chiral limit and all the new parameters that appear up to order m 2 q can be fixed 
except one; namely, the second derivative of the quark condensate with respect to 
quark masses. Some predictions of the model we are discussing involving coupling 
constants and masses of vectors and axial-vectors in the presence of masses are 

fvij Myij = fvkl My M , 

fvij M vij - fvki m vm = ~ 2 W)x ( m i + m J -m k -mi) 

flij M^ij + fij — flki Mam + fli j 
flij M Aij - f AU M Akl = ^(qq) x {rrn + rrij - m k - mi) , (22) 

where i, j, k, I are indices for the up, down and strange quark flavors. 

Three-point functions (see Figure 3) are just calculated at present in the chiral 
limit - outside the chiral limit three-point functions will be ready soon [26] . Some 




/ \ 
/ \ 
/ \ 
/ \ 



Figure 3: Infinite geometrical series which gives full three-point functions at large 
N c . The crosses that glue the vertex functions are either gy for vector or axial- 
vector sources or g s for scalar and pseudo-scalar sources. 

three-point functions have also been calculated in other large N c approaches also 
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in the chiral limit [15, 16, 17, 18] -for instance, PVV, PVA, PPV and PSP three- 
point functions, where P stands for pseudoscalar, V for vector, A for axial- vector 
and S for scalar sources. They agree fully with the ones we get in our model 
when restricted to just one hadronic state per channel. As we said before one 
could also add systematically more hadronic states per channel in our large N c 
model. 

For the moment, we have just calculated the full four-point functions -and 
the corresponding four-point vertex functions- needed for (12) and in the chiral 
limit but we are finishing their calculation outside the chiral limit too. 



3 Results 

After integrating over the four-dimensional Euclidean solid angle Qq and doing 
the limit q 2 — > as in (16), we get 

FX [Q 2 ] = 7^^ + 7^ + 7^+ ^ 



Q 2 + M 2 Q 2 + M\ Q 2 + M 2 (Q 2 + M 2 y 

+ ?a + + 1a (9-v) 

{Q 2 + M 2 A f {Q 2 + M 2 f (Q2 + M 2)3- 



This is the input needed in (15) to get B\. The explicit calculation reveals that 
diagrams with the exchange of vector and axial-vector states produce not only 
single poles as one may naively assume but two- and three-poles terms in those 
states masses too. On the opposite, scalar exchange just produce single poles 
-this is due to the different powers of momenta involved in the scalar vertices 
with respect to the spin one ones. 

The function F X [Q 2 } in (16) reproduces the large iV c -pole structure found in 
[22, 23] but including the first hadronic state in all the channels. The coefficients 
OLA, olvi OLS, Pv, Pa, lVi and 7a are known functions of M v , M\, Mj, Fq and 
the four unknowns we mentioned above. We need now to fix these inputs. 

For My, Ma and Mg, one could either identify them with the measured masses 
of the corresponding lowest lying hadronic states. That means the first multiplet 
for the vector and axial-vector multiplet while in the scalar case, the second 
multiplet is favored in view of the increasing evidence for a dynamical origin of 
the first -see [19] for a recent discussion of this issue. Another possibility is to 
use the one to one large N c relation to write the resonance masses in terms of the 
order p A CHPT couplings Lj. In the large N c limit both choices are identical. We 
prefer to take the second strategy since we want to have the order p 4 term in the 
CHPT expansion exact analytically. We use L 9 = 5.9 • 10~ 3 , L w = -4.7 • 1(T 3 , 
L 5 = 0.9 • 10~ 3 from [20, 21] which correspond to M v = 0.79 GeV, M A = 1.03 
GeV, and Ms = 1.43 GeV, which agree with the experimental masses within the 
typical l/N c uncertainty around 30 %. 
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Figure 4: We plot the function F X [Q 2 ] in three cases: the full case in (23) is the 
solid line, the OPE result in (17) is the dashed line and the CHPT result in (21) 
is the dot-dashed line. 

The four unknowns we are left with correspond to three-point functions, 
namely one for the PSS, one for the PAS and the other two for the PVA three- 
point function. These we fix as follows: we fix the one related to the PSS three- 
point function requiring that the slope at the origin is the one in (21); the constant 
related to the PAS three-point function is fixed by requiring that the residue of 
the 1/Q 2 term in the OPE in (17) vanishes. The two constants related to the 
three-point function PVA are set by requiring at the same time good matching 
between our F X [Q 2 } and the OPE result including dimension six and eight [23] 
and that this matching is between 1 and 2 GeV 2 . We checked that the depen- 
dence on the exact point of the matching is negligible. That fixes completely all 
the unknowns. 

Notice that we have analytical scale and scheme independence as was ex- 
plained in Section 2. The result for F X [Q 2 ] is plotted in Figure 4. 

We now use (15) and integrate (16) up to the matching point \x and from that 
point on we use the OPE result including dimension eight corrections [23]. This 
OPE contribution to B\ is negligible. In Figure 5, we plot the value one obtains 
for B\ integrating (16) up to the matching scale \l. Notice the nice plateau one 
gets between 1 and 2 GeV 2 . Varying all the inputs: values of L^s needed, a$, ■ ■ ■ 
between their respective uncertainties we get as preliminary result 

B\ = 0.29 ±0.15 (24) 

which is somewhat smaller but fully compatible with the one found previously 
also in the chiral limit in [12, 13, 22, 23]. 
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Figure 5: B\ plotted vs the upper limit of the integral in (16). See text for 
further explanation. 

4 Prospects and Conclusions 

The full analysis of Bk including chiral corrections is under way. The results 
of the two energy regimes explained in Section 2 are known. It is easy to get 
that F[0] = for the real case instead of F x [0] = 3 in the chiral limit. Notice 
that the value of F x [0] is a strong constraint on the value of B K . The CHPT 
calculation in the real quark masses case -F[Q 2 ] has been done to order p A [27] 
and -F[Q 2 ] remains small -below 0.15- up to energies around 0.2 GeV 2 and then 
goes negative. The other energy regime in which -F[Q 2 ] is also known is for very 
large values of Q 2 , which can be calculated using the OPE in QCD. In fact, the 
dimension six operator and Wilson coefficient are the same as the chiral limit 
one. 

One expects that higher CHPT order terms correct the behaviour of -F[<5 2 ] 
for Q 2 > 0.2 GeV 2 and F[Q 2 ] will tend to the chiral limit curve in Figure 4. The 
question is at which energy does it happens. This can only be answered using a 
hadronic large N c approximation to QCD at present. So, though we have strong 
indications that the value Bk = 3/4 has small chiral corrections as shown before, 
we have to wait till we get the full four-point Green's function in the real case to 
confirm it (12) . 

In the past few years, lattice QCD has produced many calculations using 
different fermion formulations, including preliminary unquenched studies -see 
[7, 8] for references. The recommended lattice world average result in [8] is 
Bk = 0.8llo!i3- There have also appeared chiral limit extrapolations -like for 
instance in [28] obtaining B\ = 0.32 ±0.22 or in [29] obtaining B X K = 0.34 ±0.02 
-which show a clear decreasing tendency with respect to the real case value, in 
agreement with the results found here and in [12, 13, 22, 23]. 

The QCD-Hadronic Duality result for Bk[3, 4] is very close to the chiral limit 
result above because -as already mentioned in [4]- what was calculated there, is 
the order p 2 coefficient of the chiral expansion which actually is the chiral limit 



value of B K . 
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